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ABSTRACT 

Continuing our previous QCD Hamiltonian studies in the gluonic and quark 
sectors, we describe a new renormalization procedure which generates an effec- 
tive Hamiltonian. The formulation, which is in the Coulomb gauge, provides an 
improved framework for investigating hadron structure. 



1. Introduction 

For the past few years the NCSU group has focused upon the development of 
a comprehensive relativistic, many-body description of hadron structural. This ap- 
proach has many virtues, particularly the explicit incorporation of gluonic degrees of 
freedom and a rigorous connection to QCD at short distances. Calculations^ of the 
glueball spectrum have been quite encouraging, reproducing the available lattice data. 

This talk details recent results from a new renormalization procedure for the 
Coulomb gauge QCD Hamiltonian. There are several reasons for selecting the Coulomb 
gauge. First, it is completely physical; there are no spurious gauge degrees of freedom or 
ghosts. Second, both qq and gluon-gluon interactions appear explicitly in the canonical 
Hamiltonian which facilitates contact with quark model phenomenology. Finally, there 
appears to be a natural scenario for confinement in the Coulomb gauge, as originally 
suggested by Gribova and recently formalized by Zwanzigera. 

This Hamiltonian must be supplemented with a regulator, which is applicable 
nonperturbatively and is consistent with the many-body techniques we eventually wish 
to employ. Next the renormalized Hamiltonian is constructed, which includes the cor- 
rect cut-off dependent counterterms such that physical quantities (e.g., eigenvalues) 
are cut-off independent. It is inevitable that the regulators at our disposal will vio- 
late subsets of Lorentz covariance and gauge invariance, so that the structure of the 
counterterm Hamiltonian may be expected to be quite complicated. Due to asymptotic 
freedom, however, it should be reasonable to construct these operators using pertur- 
bation theory, provided the energy cut-off is not too low. 

The result is a Hamiltonian formulation of QCD which is fully covariant and gauge 
invariant through some order in perturbation theory. To effectively apply this Hamil- 
tonian, further approximations and phenomenological considerations are necessary. An 



approach we have initiated involves making a variational Ansatz for the physical vac- 
uum state, based on the BCS "pairing" form. Quasi-particle degrees of freedom emerge 
naturally with constituent-scale masses and may be identified with constituent quarks 
and gluons. One then can construct approximate bound states of these constituents 
utilizing standard many-body techniques such as Tamm-Dancoff and RPA. 

2. Regularization and Renormalization 

Our starting point is the canonical QCD Hamiltonian in the Coulomb gauge 
V • A = 0. The dynamical degrees of freedom are the transverse gauge fields A, their 
conjugate momenta n, and the quark field if>. The canonical Hamiltonian may be 
written as 

H Ca , a = H + Hgg + Hgg + H C , (l) 

where H Q is the free Hamiltonian for quarks and transverse gluons, H gg contains the 
self-interactions between transverse gluons, and H qg contains quark-gluon interactions. 
He is the Coulomb term, 

H c = \g 2 J d 3 xd 3 y .r-y (x)/T*(x, y)J P \y) , (2) 

where J = det [V • D] is the Fadeev- Popov determinant, and the kernel K is given by 

K ab (x,y) = (x,a|(V-D)- 1 (-V 2 )(V-D)- 1 |y,6) . (3) 

The color charge density p a includes both matter and gluonic contributions, p a (x) = 
if)' (x)T a -?/>(x) + f abc A b (x) ■ II c (x), so that He represents a complicated set of instanta- 
neous quark-quark, quark-gluon, and gluon-gluon interactions. 

Perturbative schemes for renormalizing Hamiltonians typically suffer from tech- 
nical problems related to the occurrence of vanishing energy denominators. One way to 
avoiding this difficulty is the cut-off method of Glazek and WilsonH. In this approach 
one considers the basis formed by eigenstates of the free Hamiltonian Ho, with eigen- 
values E n . The theory is then regulated by suppressing matrix elements of H between 
states for which the absolute value of the free energy difference E nm = E n — E m is 
large. Specifically, we define matrix elements of the regulated Hamiltonian H(A) by 

(n\H(A)\m) = E n 5 nm + / mn (A)(n| J Ff / |m) , (4) 

where Hi = H — Ho is the interaction. Here f nm (A) is some convenient function which 
is unity for \E nm \ <C A and vanishes for \E nm \ ^> A. In Ref.i we chose 

/ nra (A) = e~ E -' K2 . (5) 

Note that, since the cut-off is defined in terms of free energy differences, it is not 
fully Lorentz covariant (though it is of course rotationally invariant). It also violates 
gauge invariance, though this issue is difficult to address directly in the present context 
since we work in a fixed gauge. The counterterm Hamiltonian will therefore also contain 
Lorentz- and gauge-noninvariant terms, which are necessary to correct for the violations 



induced by the regulator. These complications are in practice unavoidable if one wishes 
to work with Hamiltonians in a relativistic field theory. 

At this point we have a fully regulated formulation of QCD in the Coulomb 
gauge. Our next task is to remove the dependence on the cut-off parameter A by 
adding counterterms to the Hamiltonian. We thus obtain the renormalized effective 
Hamiltonian, H e s(A), which may then be analyzed nonperturbatively using many-body 
techniques. 

An elegant way of determining the counterterms involves designing a similarity 
transformation which changes the cut-off in all matrix elements of H. Because the 
transformed Hamiltonian is equivalent to the original one, physical predictions (e.g., 
eigenvalues) are unchanged. One then searches for Hamiltonians which are form in- 
variant, or "coherent," under this transformation. Such self-similar Hamiltonians will 
automatically yield cut-off independent predictions and thus represent renormalized 
Hamiltonians. Specifically, we write 

H(A') = U(A'; A)H{A)U\A'; A) , (6) 

where U is a unitary matrix constructed so that the matrix elements (m\H(A') \n) are all 
proportional to exp(— E^ m /A' 2 ). Construction of a suitable U is may be accomplished 
following the formulation of WegnerQ. We define the operator T via 

(7) 

The evolution of H under a change in A is then given by 

dH(A) 



dA- 



[T(A),tf(A)], (8) 



from which T is seen to be a generator of infinitesimal scale transformations. It is 
straightforward to show that the choice 

T(A) = [H ,H(A)] (9) 

results in the transformation properly changing the value of the cut-off_parameter in 
the regulating functions used for all matrix elements of the HamiltonianHI. 

Eqn. (|8|) can then be analyzed order-by-order in perturbation theory by expanding 
the interaction Hi = J2 P 9 P H P . At second order, for example, we find 

(n\SV 2 \m) = ^E(^ + -j=h) (e~ 2E "< E ^ A2 - e -^ m! /A' 2 ) (n|Vl(A)|/)(/|Vl(A)|m) , 

(10) 

where (n\H p (A)\m) = exp (— E 2 m /A 2 ) (n\V p (A)\m) defines the "reduced" interaction 
V p , and (n\5V2\m) = (n|T^(A')|m) — (n|T4(A)|m). This shift in the interaction explicitly 
accounts for the physics removed by changing the cut-off from A to A'. 

As discussed above, requiring form-invariance of the Hamiltonian under the sim- 
ilarity transformation determines the counterterm structure uniquely. We have imple- 
mented this constraint and determined the complete effective Hamiltonian through 



second order. The operators that appear in the counterterm Hamiltonian include, as 
expected, all possible one-body operators which are rotationally invariant (and respect, 
e.g., time-reversal invariance, charge conjugation, etc.). These are a gluon "mass" term, 
A 2 , the operators II 2 and (VjA,) 2 , and, for fermions, ipip. Each has a calculable coeffi- 
cient, except for the particular combination of II 2 and (VjA,) 2 which has the structure 
of a gluon wavefunction renormalization. For this combined operator, the variation 
of the coefficient under changes in the cut-off is calculable, and gives the gauge field 
anomalous dimension: 7^4 = — {g 2 /^Q7f 2 ) (4(7^/3), where Ca is the Casimir invariant of 
the adjoint representation. See Ref.El for full details. 

3. Work in Progress and Future Directions 

Having obtained the renormalized Hamiltonian, we can now proceed to construct 
a model for the_^ground state and analyze excitations such as glueballs, mesons and 
hybrids. In Ref£l, we focused attention on the gluonic sector and made a variational 
BCS vacuum Ansatz. This leads to a gap equation which may be solved numerically 
to give the spectral function of the quasi-particle constituent gluons. 

Work in progress includes extending the similarity analysis to higher orders, ap- 
plying this smooth cut-off treatment to the quark sector and designing improved models 
for the vacuum. At third order, one should see the running of the coupling, as well as 
the appearance of many new operators in the counterterm Hamiltonian. Going beyond 
BCS is expected to be particularly important in a combined treatment of the quark 
and gluon sectors. In addition, other composite operators, such as currents, must be 
properly renormalized which is also possible using the formalism presented in Ref.EJ. 

Finally, there are a host of phenomenological applications to be explored. Not 
only meson spectra but also form factors, electro-weak and hadronic decays, etc. can 
be predicted and confronted with available data. This comprehensive analysis should 
provide significant insight into hadron structure as well as a key linkage between fun- 
damental QCD and the successful phenomenological quark model. 
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